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We calculate the reaction cross-sections for the fusion of hydrogen and deuterium in strong 
magnetic fields as are believed to exist in the atmospheres of neutron stars. We find that in the 
presence of a strong magnetic field {B ^ 10^^ G), the reaction rates are many orders of magnitude 
higher than in the unmagnetized case. The fusion of both protons and deuterons are important over 
a neutron star's lifetime for ultrastrong magnetic fields {B ~ lO^^G). The enhancement may have 
dramatic effects on thermonuclear runaways and bursts on the surfaces of neutron stars. 
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I. INTRODUCTION: ATOMIC STRUCTURE IN AN INTENSE MAGNETIC FIELD 
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In large magnetic fields a hydrogen atom is compressed both perpendicular and parallel to the field direction. In a 
suSiciently strong magnetic field {B ^ 10^^ G), the Schrodinger equation for the dynamics of the electron separates 
into axial and perpendicular (azimuthal and radial) equations. As the potential is axisymmetric around the direction 
of the magnetic field, we expect no azimuthal dependence in the ground-state wavefunction of the electron. 

In the direction perpendicular to the magnetic field, the wavefunction can be obtained exactly [Q. This azimuthal 
wavefunction is denoted by two quantum numbers n and m. Here we take n = 0, as the n > solutions are less 
bound and therefore provide less shielding. 

The perpendicular wavefunction has the same form as the Landau wavefunction for an electron in a magnetic field: 
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A. The axial wavefunction 



Along the direction of the magnetic field, the electron experiences an effective potential. 
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Performing the integral yields 
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which for large z approaches —e^/z. The Schrodinger equation with this potential is not analytically solvable. We 
can note certain features of the desired solution. Because, Vcs is everywhere finite, both the wavefunction and its 
first derivative must be continuous. Rather than solve the equation directly, we use a variational principle, which 
constrains the ground-state wavefunction {v = 0) along the magnetic field for the given values of n and m. The index 
v counts the number of nodes in the axial wavefunction. As with the n > states, the i' > states are barely bound 
compared to the v = state. 

Looking at the radial wavefunction, we take the wavefunction along the z-axis to be a Gaussian as well: 
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We must minimize the integral, 

For this problem the integral is (using the definition of Z and m = 0), 

e^V^exp(uV2)erfc(M/V2) 
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where we have substituted u = z/an- Next, we use the definition of Z, 
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The first integral is tractable yielding the quantity to be minimized 
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with respect to a = az/au- This minimization yields a value of a^. Tab. Q lists the results for the minimization for 
several magnetic field strengths and compares them with the eigenvalues for the energy of the bound state derived by 
Ruder et al. ^ . Ruder et al. use a series of basis functions to solve the Schrodinger equation. 

TABLE I. The results of the minimization. 
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Our binding energies fall short of theirs by approximately twenty percent, because we are restricted by our trial 
wavefunction. We also tried a sum of Gaussians but this added degree of freedom did not yield significantly more 
tightly bound wavefunctions. 

Using the results of the minimization, the electron probability density is 



p(r, z) 
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where we have combined the two Gaussians in a revealing fashion. The quadrupole moment of the distribution is 
given by Q = 2a\j{a^ ~ 1)- Next we define a quantity 
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and recast the previous equation into the form 
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B. The screening potential 



When solving gravitational problems one often looks for electrostatic analogues. Here, we look for a gravitational 
analogue to an electrostatic problem. The density of the electron is constant on concentric, similar homoemoids. For 
this density distribution the potential is directly solvable M 
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In our case, we use G — — e^, ai — a2 — an, 03 — Uz and 
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Substituting these results into Eq. ^ yields 
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We change variables to simplify the integral. Using the natural units of the problem, we let r = r/an, z — z/an, and 
u = T I aj^ . This gives the new equation 
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where we again use the previous definition of a. The potential at the center of the electron cloud (r = 0, z = 0) is 
given by 
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Moving away from the origin a change of variables is useful when evaluating the integral. Let 
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The integral becomes 
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As an example we present results for B — 9.8 x 10^^ G. For this field strength an ~ 10^^^ m and ~ 10^^^ m, so 
a — 10. The range of the nuclear force is approximately 10^^^ m or 0.001a//. Fig. |l| depicts the potential in units 
of e^/a/f for this configuration. The central potential is approximately 0.25e^ /an and drops quickly in the radial 
direction. In the axial direction, the potential forms a "core" 
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FIG. 1. The left panel depicts the screening potential as a function of radius and z-position. The right panel shows the 
total potential experienced by an incoming proton. The dashed contour denotes zero potential. The other contours are 
logarithmically spaced. In the left panel the bold contour traces a potential of O.le^ /an- The contour levels increase toward 
the center with a spacing of 10^'''^'' In the right panel the bold contours trace potentials of ±10~^, 10~^ . . . e^/an- 



The total potential of the electron cloud and the proton may be approximated by the quadrupole formula 
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for large separations. 



C. The Cloud-Cloud Potential 



When we consider the interaction between the two electron clouds surrounding the protons, we must account not 
only for their electrical potential but also the antisymmetry of the mutual electron wavefunction. Because of the 
strong ambient field, we expect that both electron spins will be aligned with the field, so the spatial component of 
the wavefunction must be antisymmetric. That is. 
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with p(x) given by Eq. ^ and xq is the position of the center of the second electron cloud. 
The potential energy of the two electrons is given by {e.g. B) 
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where we have used the Gaussian form of p{x) to simphfy the expressfon for J{xo), and to obtam its approximate 
value we replace xi — X2 by xq and similarly for the other coordinates. 

To calculate the direct term of cloud-cloud potential {A{xq)) we will take advantage of the special form of the 
density distribution given in Eq. 1 1 , The direct term is in general given by 

A (xo) = J d^xip {xi - xo) $ (xi) (32) 

which is simply the convolution of the density distribution with the potential. If we perform the Fourier transform of 
the right-hand side we get 



A{xo) = J d^fcp(fc) $ (fc) e-*^ "«. 
Expressing the Poisson equation in Fourier space gives 
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Because the magnetic field induces the deformation of both electron clouds, the clouds are aligned, and they have the 
same Fourier transforms; therefore. 
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Because is a three-dimensional Gaussian, so is its Fourier transform; consequently, 
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Combining Eq. pq and Eq. p6l yields 
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Comparing this equation with Eq. p4, we get 
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Therefore, the total potential energy between two hydrogen atoms separated by x in the magnetic field is given by 



where $ (x) is simply the potential induced by the Gaussian cloud of charge (Eq. |2|) . 

Far from the atoms (r ^ aan), the interaction energy may be approximated by the quadrupole-quadrupole energy, 

Vix)^^ {a^ ~ if ^ (35 cos^ 0-30 cos^ + 3) . (41) 

where 6 is the angle relative to the symmetry axis of the atom. 

Fig. ^ and Fig. ^ depict the total potential energy between two hydrogen atoms in a magnetic field for the same 
magnetic-field strength as Fig. |^ (i? = 9.8 x lO-'^^G). A comparison of the two figures illustrates that the exchange 
term provides a slight attractive force between the two electron clouds, because of the anticorrelation of the clouds. 
At large separations, both potentials are well approximated by the quadrupole-quadrupole formula (Eq. ^T|) . 
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FIG. 2. The figures depict the total potential energy between two magnetized hydrogen atoms excluding the antisymmetriza- 
tion energy. For the left panel, the contour spacing is the same as in right panel of Fig. |^ The right panel illustrates the 
potential along the axis of the magnetic field. 
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FIG. 3. The figures depict the total potential energy between two magnetized hydrogen atoms including the antisymmetriza- 
tion energy. For the left panel, the contour spacing is the same as in right panel of Fig. |l|. The right panel illustrates the 
potential along the axis of the magnetic field. 



II. ESTIMATING REACTION RATES 



In a fluid state, there will be three possible reaction channels, 

• proton- proton dominates in hot, totally ionized gas 

• proton-atom dominates in nearly completely ionized gas 

• atom-atom dominates in neutral and partially ionized gas 

For the flrst channel, we can use the standard thermonuclear reaction rates {e.g. jsj). For the latter two channels 
we must include the screening potentials that we calculated in the previous section to determine the potential wall 
through which the interacting particles must penetrate. 



A. The transmission probability 



In the WKB approximation, the probability to traverse through a potential wall is 
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where B12 is the magnetic-field strength in units of 10^^ G, u is the dimensionless radius r/an and £ and V are the 
dimensionless energy Ean/e'^ and potential. 

For the proton-atom channel, the potential includes both that of the nucleus V = 1/u and the surrounding electron 
cloud (Eq. HI). At large distances from the nucleus, u » a, the total potential is well approximated by the 
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quadrupole (Eq. [2^ ). For the atom- atom channel, the total potential includes contributions from the proton-proton, 
proton-electron and electron-electron potentials (Eq. which is well approximated by the quadrupole-quadrupole 
formula (Eq. ^l]) for large separations. 

To calculate the transition probability, we use these quadrupole formulae to approximate the potential for u > Aa 
and for u < 1/2, we approximate the potential energy between the electron clouds and the electron clouds and the 
protons by their central values. This both speeds the calculation and reduces the numerical error. 

Fig. |4| traces the transmission probability for protons to interact with atoms and atoms to interact with atoms at 
zero relative energy as a function of angle and magnetic field. In the atom-proton case, the protons can most easily 
penetrate through the mutual potential barrier along the axis of the magnetic field and the penetration probability 
increases markedly with the strength of the magnetic field. In the atom-atom case, we see that the maximum 
transmission probability occurs at an angle to the field direction and that with antisymmetrization of the electron 
density the transmission probability increases dramatically. For the reaction rate estimates that follow we will account 
for the antisymmetrization energy of the two electron clouds. 




FIG. 4. The left panel depicts the transmission probability as a function of angle and magnetic field for a proton and an 
atom to interact at zero relative energy. The right panel depicts the same probability for two atoms. The solid lines trace 
the probability if the antisymmetrization energy of the electrons is considered. The dashed lines show the probability without 
antisymmetrization. 

To translate this transmission probability into a cross section, we must average |Tp over a sphere and include the 
appropriate S-factor for the reaction where S{E) is defined as 

S{E) ^ aE\Tf ^ Sail + SiE) (45) 

In this way, the strong energy dependence of the reaction cross section is removed. For the reaction ^II(p, e+i/)D, 
So = 4.38 X 10-25 MeV-barn and Si = 11.2 MeV"^ at low energies m . The reaction of the less abundant deuterons 
with protons has a much larger S-factor of Sq = 2.5 x 10"'^ MeV-barn and 5*1 = 27.8 MeV"^ 

Given these definitions, the reaction proceeds at a rate of 

ri2 = {1 + 6i2)~^nin2 (av) (46) 
/ o \ 1/2 

= (l + ^i2)-'nin2 — 5o(l + 5ifcT)(fcr)-i/2|j^|2 (47) 
\fnTj 

where fi is the reduced mass of the reactants, and rii, n2 are their number densities. ri2 has the units of reactions per 
unit time per unit volume, so we can define a typical timescale for a reactant to be consumed 

Ti2 = ni/ri2. (48) 
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We will use this timescale to assess the effectiveness of the screening in catalyzing the nuclear fusion reactions. We 
also account for the increasing excitation of the gas as the temperature increases and the onset of thermonuclear 
reactions above several million degrees. 



B. The ground-state fraction 



The screening is much less effective if the electron is in an excited state, so we estimate the fraction of atoms in 
the ground state by first calculating the ionization equilibrium and then the fraction of neutral atoms in the ground 
state. 

Lai and Salpeter give the form of the Saha equation for hydrogen atoms, electrons and protons in equilibrium 
in the presence of a quantizing magnetic field. Throughout this formalism, we use the natural units of the problem, 
i.e. T is the temperature in units of 3.15 x 10^ K, M is the mass of the system in units of the electron's mass (1840 
for hydrogen and 3670 for deuterium), b is the strength of the magnetic field in units of 2.35 x 10^ G and Ug is the 
number density of the gas in units of 6.76 x lO^^cm^'^. 

We first look at the unexcited hydrogen atom. For the partition function of the ground state in a quantizing 
magnetic field, Lai and Salpeter 0] give 
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where E{H) = —0.16/^ (the ground-state energy of the atom). 
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Here we have explicitly integrated to iiTj^max, so we replace M'^ of Lai and Salpeter Q with A/", 
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and M'^ is as given by Lai and Salpeter , 
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As /4r_Lmax ^ cxD, M" M'^ and we recover the Lai and Salpeter result, if^max is the upper limit on the 
perpendicular momentum for the given state. The electron clouds of neighboring atoms should not overlap; otherwise, 



the gas would become pressure ionized. Therefore, we take the size of the state, Rk = K±/b < Rg = Ug 
defining condition on ii'_Lmax- We obtain 
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The total partition function of the neutral atom is given by 
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where and Zm are the partition functions for excitations of the v and m quantum numbers respectively. Lai and 
Salpeter Q argue that the z^{H) ~ 1 as these states are hardly occupied relative to the ionized, m > and ground 
states. For the contribution of the m > states to the partition function, they obtain 
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where we have several additional auxiliary definitions: 
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and M'^^^ is as given by Lai and Salpeter [Q, 
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M±m is given by the following relation 
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and we use the following additional definition 

T„ ~ 0.64/,„(Mj_„, 
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The ratio of the number of atoms in the ground state to the number of neutral atoms is given by 

Aground (-^) .^ground (-^) 1 
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Next we calculate the ionization-recombination equilibrium. 

X{H) Z{H) 



Zra{H) ' 



XpXe 



Z{p)Z{e) 



2tt 
X exp 



1/2 



tanh 



V 2MT 



\Em 

T 



Zrn{H), 



(58) 



(59) 



(60) 



(61) 



(62) 



(63) 



(64) 



(65) 



(66) 



where X{H) — n{H)/ng, Xp np/ug, X^ — ne/ug are the number density fractions of the different species. 

Combining Eq. p3 and Eq.^q yields the fraction of "shielded" nuclei as a a function of temperature, density and 
magnetic-field strength. Fig. ^[depicts the fraction of unexcited hydrogen atoms in the gas as function of temperature 
for several field strengths and two densities. 
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FIG. 5. The ground-state fraction as a function of temperature, density and magnetic field. The left panel shows the neutral 
fraction as a dashed line and the unexcited fraction as a solid line for p ~ 1 e cm" and B = 10^^ lO" and lO^'' G. The right 
panel is for a density p ~ 1000 g cm~^. 



C. Thermonuclear reactions 



We parameterize the thermonuclear reaction rates {e.g. Q) by 

rpp = 3.06 X 10-3^cm3sec"^n2rg"2/^exp(-33.7irg"^/^) 
rpD = 3.28 X 10-^^cm^sec-^npnDTQ^^^exp{-37.nTQ^^^). 

The timescale for the exhaustion of a particular reactant becomes 



^thcrmo ~r ^magneto 



(67) 
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Fig. 1^ shows the reaction timescale for the consumption of hydrogen and deuterium in the reactions p(p, e^v)D and 
£'(p, 7)'^He respectively for a magnetic field of lO^^G. 
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FIG. 6. The two panels depict the reaction timescale for the reactions p{p,e^v)D and D(p, 7)^He for B — lO^^G over a 
range of temperatures and densities. The dashed contour traces r of one second. The solid contours trace locii of timescales 
ranging from one year to 10^^° years with a factor of ten in between each contour. 



Even in this very strong magnetic field, the p-p reaction proceeds only very slowly below temperatures of one million 
degrees; however, over millions of years, the hydrogen gas would be processed to deuterium and then to helium in such 
a strong magnetic field. It would provide a steady source of energy, while eroding the storehouse of hydrogen which 
could potentially fuel a thermonuclear runaway. Relatively, the second reaction proceeds instantly with timescales of 
less than one year for the interesting range of densities and temperatures. 

For the weak fields depicted in Fig. M only the deuterium reaction proceeds at a significant rate. 



D(p,7)-'He B^lxIQi* G D(p,7)-'He B^lxIQi^ G 




101 102 10^ 10-* 105 10^ 10^ 108 10^ 10^0 101 102 10^ 10-* 105 106 10^ 108 103 IQIO 

p(g cm"-) p(g cm"-) 

FIG. 7. The two panes depict the reaction timescale for the reaction D{p,'y)^ He for B = lO^^G and 10^'^G over a range of 
temperatures and densities. The dashed contour traces r of one second. The solid contours trace locii of timescales ranging 
from one year to 10^" years with a factor of ten in between each contour. 
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III. DISCUSSION 



We find that in strong magnetic fields {B ^ 10^^ G), the cross-section for nuclear fusion is dramatically larger than 
in the unmagnetized case. For these strong fields, deuterons fuse to ■^He over short timescales (<; 10^ yr) for the 
density and temperatures expected on the surface of a neutron star. Because of the inherent weakness of the p — p 
interaction, the fusion of protons to deuterium is only important over cosmological timescales for ultrastrong fields 
{B ^ 10^^ G) in spite of the large enhancement in the cross section of this reaction. 

For larger atoms {Z > 1), we expect that reaction cross-sections will also be larger in the presence of an intense 
magnetic field. However, the shielding is unlikely to be as effective as for the Z — 1 case, because additional electrons 
must occupy m > levels which are much less effective at screening the nuclear charge. 
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